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PACS 03 . 65 . Yz - Decoherence; open systems; quantum statistical methods 
PACS 03 . 65 . Fd - Algebraic methods 

Abstract -For a certain class of open quantum systems there exists a dynamical symmetry 
which connects different time-evolved density matrices. We show how to use this symmetry for 
dynamics in the Liouville space with time-dependent parameters. This allows us to introduce a 
concept of generalized coherent states in the Liouville space (i.e. for density matrices). Dynamics 
of this class of density matrices is characterized by robustness with respect to any time-dependent 
perturbations of the couplings. We study their dynamical context while focusing on common 
physical situations corresponding to compact and non-compact symmetries. 



Introduction. — The concept of coherent states plays 
a very important role in quantum physics. Introduced 
by Schrodinger for quantum harmonic oscillator [T| , it re- 
ceived further application with the birth of quantum op- 
tics in the 60's [2j[3]. Further developments of the concept 
of coherent states are associated with the generalized co- 
herent states defined for any Lie algebra, with the most 
important examples given by the su(2) and su(l, 1) alge- 
bras. Coherent states have extremely wide applications in 
physics and mathematics, reviewed e.g. in [4j|6]. 

There are several definitions of the coherent states. The 
oscillator coherent state can be defined either as (i) a state 
which minimizes an uncertainty relation, or (ii) an eigen- 
state of an annihilation operator, or (iii) a state obtained 
from the vacuum state by the action of the displacement 
operator. These three definitions are equivalent for the 
harmonic oscillators while they are not for the more gen- 
eral (generalized) coherent states associated to some non- 
abelian algebras, other than a Heisenberg algebra. The 
construction of coherent states associated to a Lie alge- 
bra includes three ingredients: (i) an algebra g with the 
representation space; (ii) a vacuum state of ladder oper- 
ators in this representation space. This state is supposed 
to have an invariant subgroup defined by some subalge- 
bra h € g\ (iii) a concrete representation of a group G. 
It follows from this that generalized coherent states have 
profound geometrical interpretation: they are labeled by 
the points of a homogeneous (coset) space G/H which in 
many important physical situations has a structure of a 
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Kahler manifold. 

Those quantum-mechanical systems which have a dy- 
namical symmetry given by G (e.g. the generators of G 
commute with the Schrodinger operator id t — H, and not 
just with the Hamiltonian) can be described by a classical 
dynamical system on a coset G/H. One of the most im- 
portant properties of the generalized coherent states exists 
due to their Lie group structure: if the initial state of a 
quantum system is a (generalized) coherent state, it will 
remain so at any subsequent evolution time (see eg. [7]). 
This key property will be important for the construction 
below. 

The developments briefly outlined above concern the sit- 
uation of isolated quantum systems. In this letter we pro- 
pose a construction scheme of generalized coherent states 
for a certain class of quantum systems with dissipation. 
The most convenient formulation is in terms of the Liou- 
ville space on which the dynamics of the reduced density 
matrix is defined by the Liouville operator C The dynam- 
ical symmetry in this case is associated with a symmetry 
of the density matrix, so that there is an algebra G of op- 
erators commuting with d t — C. We thus define generalized 
coherent states for an open system as Liouville coherent 
states (LCS). These states allow for a new interpretation 
of the Liouville dynamics and simplify the calculation of 
typical quantities of interest. In this paper we first intro- 
duce the very concept of LCS and then we demonstrate 
it on several important examples: the simplest non-trivial 
ones when the algebra G is isomorphic to su(2) algebra, 
which happens for spin-boson-type models, or su(l, 1) al- 
gebra, which corresponds to models of harmonic oscillator 
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type. A defining property of these density matrices is a 
robustness of evolution with respect to any time-pendent 
driving. 

Dynamical symmetry approach to open quan- 
tum systems. — In general, the formal solution of the 
evolution equation for an open system is given by the time- 
ordered exponent p(t) — Texp(Ci(t))p(Q), where £7 is 
a Liouville superoperator (that is, it acts on operators) 
evaluated in the interaction picture. When focusing on a 
subsystem, the reduced density matrix is evolved in time 
by the Lindblad generators which are constructed using 
the eigenoperators of the subsystem according to the usual 
technique (see e.g. (8|9]). It is possible that these superop- 
erators fulfill some Lie-(super)-algebra identities [lOj. This 
is the starting point of our construction. We demonstrate 
that a large class of physically relevant systems indeed 
satisfies this assumption. 

We consider an evolution of the reduced density matrix 
in the Lindblad form 



simplest case of the su(2) and su(l, 1) algebraic structure, 
defined by the commutation relations [£_,L+] = 2a Lq, 
[Lo,L±] — ±L±, where a — ±1 refer to the sit(l, 1) and 
su(2) cases, the function /+ can be shown to satisfy the 
following Riccati-type equation {j G {+, z, — }), 



/+ - 7«(*)/+ - <n-i*)f+ - 7+(*) = 0, 



(4) 



while the remaining functions are determined as 
/*(*) = Jj(7*(*) + 2o7_(t)/+(t))d* and /_(*) = 
J 7_(t) exp(f z (t))dt. This equation can be solved (ei- 
ther analytically or numerically) for any functional form 
of 7± (t) , 7 2 (t) , in particular for the case when they do not 
depend on time. The solution in such a case reads 



/±W ~ cosh(i^)- 
f a (t) = [cosh(tD) 



(j±/D)swh(W) 



(7o/2£>) sinh(iD) ' 
- ( la /2D) smh(tD)}- 2 , 



(5) 



C(t)p(t) 



- h \H{t),p{t)] 



(1) 



J2^(WA jP (t)A] - A)A, jP (t) p{t)A]A 3 ), 



which is the only form compatible with the positivity of 
the density matrix. Here Aj are eigenoperators of the 
subsystem. In general, j(t) = f Q drRe[Y(t — r)], where 
Y{x) — / °° dujuj^ 1 p(oj) exp(-iujx) is related to the spec- 
tral density p(ui) which contains information about sub- 
system-reservoir coupling. Introducing the following no- 
tation: (A <Ei B)p := ApB, one defines the action of the 
individual terms in on the density matrix as Lip (e.g. 
Lip = ApA\ L 2 p = A^Ap, L 3 p = pAA\ etc.). 

Now we assume that the components of the Lindblad 
operators obey Lie-algebra type relations 



[Li, Lj]p — fijkLkP, 



(2) 



where fijk are some numbers (structure coefficients) . Our 
approach is also valid if fijk are proportional to some op- 
erator which commutes with all the other operators. The 
assumption ^ is satisfied for a surprisingly large num- 
ber of systems, basically for all the standard Liouvillians 
discussed in literature, and it was shown to be a rather 
general property of Lindblad-type dynamics 10 . 

Provided that the relations ^ are satisfied, the time- 
ordered exponent for any time-dependent coefficients jj (t) 
is an element of a Lie group corresponding to the Lie alge- 
bra of superoperators. Using the disentangling technique 
(see e.g. [9]), the time-ordered exponential can be trans- 
formed into a product of ordinary exponentials 



p(t)=l[cxp(f J (t)L J )p(0)- 



(3) 



The relation between the functions fj(t) and the func- 
tions Jj(t) can be derived easily. In particular, for the 



where D = ((70/2) 2 - crj+j^) 1 / 2 . 

The disentangled form (J3J) allows a direct evaluation of 
any dynamical quantity of interest, such as the entropy 
S = — Tr(plogp) or the purity of the density matrix 
Tr(p 2 ), for arbitrary time-dependent couplings. 

The use of the Lie group structure greatly simplifies 
computations of any correlation functions. For example, 
the quantity Tr(Ap(t)) with A £ g can be reformulated 
as d\Tr(exp(XA) p(t))\\ = o; the expression in the brackets 
can be evaluated using the group multiplication rules JHJ. 

Liouville coherent states (LCS). — We proceed 
with the Liouville space formulation. A Liouville space C 
is defined as a direct product of two Hilbcrt spaces, C = 
H^T-L, corresponding to the left and right vectors in terms 
of superoperator notations. Thus, a vector \A) G C corre- 
sponds to an operator as follows \A) = ^2A mn \m,n) 
Em n Ann l m )( n li where the sum runs over orthonormal 
bases \m) and |n) of the spaces T-L and % of dimensions 
d, d. To define a geometry on the Lioville space we use the 
usual scalar product defined by (A, B) := Ti[A^B] (cf. [9]). 
This scalar product allows us to introduce the bra-vectors 
for each ket-vector in the Liouville space. The equation 
for the density matrix p = —i[H, p] gets converted in this 
notation into \p) = —iH \p) where H = H — H . This map 
is known as a Choi-Jamiolkowski map between states and 
operators and it may be used to extend certain concepts 
known for states in the Hilbert space to states in the Liou- 
ville space, which correspond to operators in the Hilbcrt 
space: the operators A = k aj ^ \ j)(k\ are mapped to 
states \A) = J2j,k a J,k \^ k )- 

Definition of LCS: We define a Liouville coherent 
state as a generalized coherent state in the Liouville space. 
The operators creating generalized coherent states thus 
act between two Hilbert spaces constituting the Liouville 
space. The general condition that the initial density ma- 
trix be in the class of coherent states means that some 
linear combination of generators of the Lie algebra anni- 
hilates the state J7 ^Li\p) = and thus defines a station- 
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ary subgroup H. We explicitly demonstrate this for the 
SU{2) and SU(1, 1) cases, in which H = U(l). 

Geometry of dissipative dynamics. Since gen- 
eralized coherent states are labeled by the points of the 
cosets G/H, where if is a stationary subgroup, the ge- 
ometry of this coset space has profound influence on the 
dynamics. The main property of this type of dynamical 
evolution can be formulated as follows: if the Liouville op- 
erator is an element of the Lie algebra of superoperators, 
and if the initial density matrix is a coherent state (in 
the sense that the corresponding vector in the Liouville 
space is a generalized coherent state) , irrespectively of the 
precise form of the time-dependent coefficients 7j(i) the 
subsequent dynamics will remain in this subclass of the 
states (density matrices). This is a direct consequence of 
the group multiplication property. The overlap between 
two such LCS is given by the cigcnfunctions of the Laplacc- 
Bcltrami operator on G/H |5;|. The distance function on 
this space is a measure of Lodschmidt echo (TT1 12 . The 



area 2- form on G/H corresponds to the imaginary part of 
the metric and describes the geometric (Berry) curvature 
of the quantum evolution. Generally, the topology of the 
coset G/H can be nontrivial. 

Examples. — We consider two types of examples 
which have great physical importance. One example refers 
to the case when the Liouville operator is a combination 
of generators of a compact group, SU(2), while the other 
example is a realization of a non-compact, SU(1, 1) group. 
For both su(2) and su(l, 1) algebras the coherent state is 
built upon a basis state (which we denote by |0)) using a 
shift-operator D(Q, such that \Q := D(Q |0), and 



10 := 



10) , (6) 



where rj = —a log (l — er|C| 2 ) and the relation between 
i and ( reads C = tanh |£|e iv (with £ = |£|e iv ) for su(l, 1) 
and C = - tan(9/2)e-^ (with £ = -(Q/2)e~ lip ) for su{2). 
Modifying the arguments of [5] for the present case, we 
may write down the resolution of identity on the Liouville 
space in the following form: 



1 = 53/ d »k(0\k,o{k,c\ 



(7) 



where k runs over all representations of the underlying 
dynamical symmetry group and p,k(0 is the group mea- 
sure [5]. 

The action of the evolution superoperator on the coher- 
ent state can be easily calculated starting from its disen- 
tangled form ([3]). We introduce the operator 

O :=e f + L +e foLo e f - L -e^ L +- CL -, (8) 

which can be rewritten in a product form 

O = e g+L+ e gaL 0e g-L- 

= D(g+) ■ e <*'-'»> £ » ■ exp([ 5 _ + g* + e^\ L_), (9) 



in which the functions gi are governed by equations 

.9+ = /+ + Ce /o A~\ e 90 = (1 - <r\C\ 2 y a e fo A~ 2 , 

5 _=-C* + (l-a|C| 2 )- CT /-A- 1 , (10) 

where A = 1 — er/_£. As expected, the evolution of an 
initial coherent state \Q leads to a coherent state \g+), 
described by an equation 



Texp£(t) |0 = c(t,g ,g+) \g+) 



(11) 



with the prefactor c(t,go,g+) = c(t) ■ e( 9o ~ I? ' fe , where k is 
a representation-dependent constant, defined by Lq |0) = 
k |0), and the factor c(i) is determined solely by the dy- 
namics within the U(l) sector. Note that the relation be- 
tween <? + and C is a Mobius transform, under which circles 
in the complex plane of Q are mapped onto circles in the 
plane of g + : a circle £ = \(\e tv is mapped to a circle of 
radius R and origin z given by 



R 



s/0| 



i-|/-l 2 ICI 5 



ofo . 



ICI 2 /- 



i-|/-l 2 ICI 5 



(12) 



Consideration of a mapping of circles is a convenient way 
to think about the evolution of coherent state parame- 
ters. Namely, one can see that while the radius of a circle 
remains approximately proportional to the radius of the 
initial circle, the origin shifts to some finite value; there- 
fore the "average" value of a coherent state parameter at 
a later time is shifted from zero. 

Spin-boson-type models: SU(2) LCS. For the 
case of spin-boson-like models in the RWA the system is 
described by a Lindblad equation of the form ([TJ , in which 
we set H(t) = a z ■ fi(t)/2, A + = a + , A_ = cr_, 7+ (f) = 
j(t)n and 7_(t) = j(t)(n + l), with n = (exp(hQ/k B T) - 
1) _1 and y(t) being system-specific. Defining 

L + p = (7 + /J(7~, L~ p = ct~ p<r + , 

L°p=±(a z p + pa z ), R P = l -{o*p~po z ), (13) 

and using the fact that cr ± (T =F = 4(1 ± cr z ), we obtain the 
following algebra: 



[L+,L-}=2L°, [L",L ± ]=±L ± 
[L^i?] = [L°,R] = 0, 



(14) 



which has a structure of the direct product su(2) x u(l) 
and consequently the supergroup operator can be disen- 
tangled using the procedure outlined above. Therefore 
the formal operator solution of the evolution equation 
p(t) = Cp(t), with 



£ 



-ifl(t)R + ~{(t)nL + + 7(t)(n + 1)L~ 



7(«)i°- S 7(*)(2fi + l), (15) 
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reads 



Pit) 



e f+L + e foL° e f-L- e - 



iR f ndt e -(n+l/2) / jdt 



P(0), 



(16) 

where f±, z {t) satisfy the Riccati equation Q. 

Any matrix, namely the density matrix can be decom- 
posed into irreducible parts classified by the Casimir op- 
erator L 2 = \ (L+L~ +L~L + ) + L Z L Z , which can be 

explicitly expressed as L 2 = | (p + cr z pa z ). In fact, us- 
ing the basis of 2 x 2 matrices cto = 1, cr±, <t z , we can 
write p = p J= Q + p 3= i/2, where p J=Q = p+a + + p-<J- 
and Pj=i/2 = PfO- + a z) + PiO- — °z)- The most gen- 
eral density matrix p for a su(2) system with spin 1/2 
can be written as a linear combination of spin coherent 
states p = c |0, +) + c* |0, -) + ci(C) |l/2; C> , C > 0, 
where |0,±) corresponds to a± (the subspace with j = 0) 
and |1/2;C) is a coherent state with parameter £ built 
upon the basis state \j = 1/2, m = —1/2), corresponding 
to (1 — a z )/2. The evolution operator acts indepen- 
dently on the two subspaces. Obviously, the coherent- 
state structure on the subspace j = is irrelevant, the 
only nontrivial evolution there being the evolution of the 
coefficients c. From the condition Tr p = 1 we derive 

<=i(0 = v / i + KF/(i + C). 



For the sake of a particular example, we follow 13 and 



consider a two- level system with level splitting u>, inter- 
acting with a one-mode bath of frequency k = uj + A. 
The non-Markovian evolution can be treated by a Lind- 
blad equation, however with time-dependent terms. A 
time-dependent energy splitting Q(t) — ^X(t) and a time- 
dependent dissipation j(t) = 5RA"(t) are given by 



X{t) 



iA't 



iA't 



where A' = ^A 2 +g 2 , g' = A' - A, s = g'/{A + A') and 
r = [u + (A + A')/2]/[u) + (A - A')/2]; we also set n = 
0. As 7+ is absent, we have /+ = 0, and the solution 



a coherent state can be deduced from Eq. (11 1 



Thanks to the absence of dissipation (the bath has one 
mode only), f—(t) is a periodic function of time and there- 
fore the coherent state parameter is periodic as well. As 
mentioned above, a circle in the plane of the initial co- 
herent state parameter £(0) is mapped onto a circle in the 
plane of coherent state parameter £(i) at time t. The "pul- 
sation" of the parameters R and z of Eq. (12 1 is plotted 



in Figure [T] while the time evolution of the purity of the 
density matrix is shown in Fig. [2j 

Harmonic oscillator-type models: SU(1,1) LCS. 

— In analogy with the LCS for spin-boson type mod- 
els, we can study models of a harmonic oscillator type, 
which lead to a non-compact symmetry group SU(1,1). 
A generic model of a harmonic oscillator in contact with 
a reservoir is described by the Lindblad equation of the 
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Figure 1: Time dependence of the parameters R and z, intro- 



duced in Eq. ( 12 l , for the su(2)-type model defined by Eq. (171 
(with parameters cj — 2 and g — 1). The initial coherent state 
parameter ( lies on a circle \(\ = 1/2. The evolution is plot- 
ted for various values of detuning A. Both R(t) and z(t) are 
periodic functions of time. 

form Q (c.f. 9 ]) with H(t) = u(t)ata, A + = a\ A_ = a, 
7+ = 72 (t)n and 7_ = 7i(t)(n + 1), where n is the equi- 
librium thermal occupancy of the oscillator and couplings 
7i,2(i) are system-specific. To proceed we denote h = a) a 
and introduce the following superoperators: 

K-p = apa i . K + p = a)pa 

K oP=l,{np + pn + p), B.p = hp~ph, (18) 



(17) which satisfy the su(l, 1) commutation relations 



[K-,K+] 



2K , [K ,K±] 
[R, K a ] = 0. 



±K± 



(19) 



to the corresponding Riccati equations simplifies greatly: K—K. 
f_(t) = -f (t) = f* j(t)dt := r(t). The time evolution of 



The Casimir invariant is given by [Kq 



l/2(K + K_ + 

4- )\p = \{—p + n 2 p + pn 2 — 2hpn). The algebraic 
structure in this case is su(l, 1) x u(l) and we may again 
disantagle the evolution operator using the method pre- 
sented above. The evolution equation p(t) — Cp(t) with 

C = [7i(t)(ft + f) - 72(t)n] - iu(t)R 

+ 2 7 i(i)(n + l)K- + 2 l2 (t)nK+ 

-2[ 7 i(t)(n + l)+ -y 2 (t)n]K (20) 

is solved by 

p(t) = e f+K+ e foKo e f-K„ 

e -»-R/w<ii e /[7i(t)(fi+l)-72(t)ft]dtp/Q^ (21) 

where f± satisfy Eq. Q for a = 1. 

The space of density matrices can be decomposed into 
a direct sum of irreducible subspaces of the underly- 
ing su(l, 1) algebra. Namely, the operators of the form 
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basis of these coherent states as 



1.0 



Figure 2: The time evolution of the purity Trp 2 (t) of the 
density matrix p for a su(2)-type system (upper panel) defined 
by Eq. (17 1, and a su(l, l)-type system (lower panel). The 



initial state of the evolution is a LCS with initial parameter £ 
(resp. £o) with different initial purities. The su(l, 1) system 
is characterized by time-dependent decay rates 71,2(4) = 7(1 + 
cos(87t))/2. Full lines describe the evolution when c m = in 
Eq. ( |22[ ) for m > (diagonal matrix). Dashed lines involve a 
small admixture of a coherent state |l,£i). 



|n + m)(n| (or \ri)(n + m\) belong to an irreducible sub- 
space, according to the action of the Casimir opera- 
tor on them K 2 \n + to) (n\ = k(k — 1) \n + m)(n\ with 
k = h (1 + to) (similarly for \n)(n + m\). Thus the space 
of density matrices is decomposed into a direct sum of 
subspaces, corresponding to discrete-series representations 
of the su(l, 1) algebra. We define \m;n) as \m + n)(n\ 
for to > and \n)(n — m\ for to < 0; then K a \m;n) = 
(k + n) \m;n), with n > 0. The evolution operator acts 
independently on each of these subspaces. 

The coherent states are introduced separately for each 
irreducible subspace m according to the recipe presented 
above in this paper. Thus we parametrize each subspace to 
by coherent states |to; £ m ) (£ m being the coherent state pa- 
rameter) built upon the basis state |to; 0). A general den- 
sity matrix can be expressed in terms of the over-complete 



m>0^ 



(22) 



where c m (C) are appropriate time-dependent functions 
and dp, m is the group measure (see Eq. Q). The co- 
herent state parameters £ m are all time-dependent; their 
evolution is given by Eq. (111. In particular, we can con- 



sider a density matrix being a purely coherent state within 
each m-sector. There are, however, constraints on the ad- 
missible set of parameters c m ,Q m , in order for the linear 
combination to indeed represent a physical density matrix. 
The simplest one reads 1 = Trp = (1 - |Co| 2 ) _1/2 [c /(l - 
£0) + c.c.]. Other constraints, such as < Tr p 2 < 1, 
are more involved; as a result not every purely coherent 
state is physically allowed. Interestingly, the decomposi- 
tion of an arbitrary physical density matrix may contain 
even these unphysical states. Because of this feature no 
simple parametrization of a coherent state in terms of its 
purity exists, contrary to the su(2) case. 

In order to demonstrate the LCS technique in the 
su(l, 1) case on a particular example, we focus on a har- 
monic oscillator interacting with a bath in such a way that 
the decay rates 71,2 (0 are both equal to a periodic func- 
tion 

7iW =72(t) =7(*) =7(a + cosrt), (23) 

where a is a free parameter. The time evolution of the 
parameters R and z of Eq. ( 12 ) can be seen in Fig. [3| 
while the evolution of the purity Tr p(t) 2 is captured in 
Fig.[2j 

Conclusion. — In this paper we introduced the no- 
tion of Liouville coherent states, an analogue of general- 
ized coherent states for the density matrices. This concept 
identifies a class of density matrices whose time evolution 
is robust with respect to any time-dependent driving: if 
the initial density matrix belongs to this class it will re- 
main so for any t > 0. We demonstrated this on several 
physical examples, involving compact and non-compact 
dynamical symmetry groups. Many geometric and topo- 
logical properties of these states deserve further study. As 
these states form an overcomplete basis, they can be used 
as a platform for investigating a more complicated Liou- 
ville dynamics analogously to the way in which the usual 
coherent states arc used in path-integral formulation of 
the quantum dynamics. 
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